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The equal-spin pairing (ESP) state, the so-called A-like phase, of superfiuid 3 He in aerogels is 
studied theoretically in the Ginzburg-Landau (GL) region by examining thermodynamics, and the 
resulting equilibrium phase diagram is mapped out. We find that the ABM pairing state with 
presumably quasi long-ranged superfiuid order has a lower free energy than the planar and "robust" 
states and is the best candidate of the A-like phase with a strange lowering of the polycritical point 
(PCP) observed experimentally. 
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PACS numbers: 

It is an important subject on liquid 3 He in aerogel 0,0 
to understand whether or not the pairing symmetry and 
other details of superfiuid states are changed due to im- 
purity effects provided by the aerogel background. At 
present, it will be clear that the B-like phase is, just 
as in bulk 3 He, the BW pairing state with superfiuid 
long range order (LRO), while it is a matter under de- 
bate whether or not the equal-spin pairing (ESP) state, 
the so-called A-like phase, near the superfiuid transition 
temperature T C (P) is the ABM state. 

By assuming that the orientation of the 1- vector in the 
ABM pairing state should lose its LRO due to the impu- 
rity disorder, Volovik has speculated that the A-like 
phase in zero magnetic field should be a glass state with 
the ABM pairing but with no superfiuid rigidity (i.e., 
p s (q = 0) =0). Fomin Q has argued with no free en- 
ergy calculation that another ESP state, "robust" state, 
should be the A-like phase in aerogel. Among experimen- 
tal suggestions on the A-like state, the lowering 0, 0, Q 
of the polycritical pressure (PCP) at which T C (P) and 
the A-B transition line Tab(P) merge is a key observa- 
tion: Although the range T c — Tab of the equilibrium 
A-like state at a fixed pressure P above the bulk PCP 
is narrower in aerogels with lower porosity, the pressure 
range over which this state is detected seems to be wider 
with decreasing the porosity 0, Q- In this paper, we ex- 
amine candidates of the A-like phase in equilibrium by 
calculating the superfiuid free energy. We find that, due 
primarily to effects of quenched disorder on the ampli- 
tude of the pair-field (superfiuid order parameter) A^i, 
the ABM state is the best candidate of the A-like state 
with the above-mentioned properties. 

To derive the Ginzburg-Landau (GL) free energy func- 
tional microscopically, we start from the quasiparticle 
hamiltonian with an effective interaction mediated by 
spin fluctuations and an impurity term 7ii mp = 
J d 3 r J2 a V'JC 1 ') u ( r ) V'cr(r), where u(r) is an impurity po- 
tential with zero mean. For a moment, impurity scat- 
tering effects inducing inhomogenuities of A^i will be 
neglected. If the impurity effects on the quasiparticles 
are included just through their relaxation time, the GL 
functional is, as usual, expressed by assuming a P -wave 



pairing interaction as the sum of Tio and 7i gr , where 



7~Lo[A p ,k] 



and 



{ocA^A^ + fclA^A^l 2 + ^(A^i) 3 

PsA^AyjAnjAvj + PaAhjAvjAjjjAhj 
^A^A^iA^jA^j), (f) 



H gr {A p , k }= A 



-KtdidjA^-^-^A^ 



(2) 



Here, J r = fd 3 r, a = A(0)[ln(T/T c0 ) + ^(1/2 + 
(47TTT)- 1 ) - ^(1/2) ]/3, K\ = K 2 = 2iV(0)£g/5, and 
Zo = vf v / -^ (2) (l/2+ (4^Tr)- 1 )/6/(4^T). Further, T c0 
is the mean field transition temperature in bulk at a fixed 
P, and r _1 = 2-7rA(0)( |zt p _ p /| 2 ) ang is a relaxation rate, 
where the overbar and (• • -) a ng denote the averages over 
the randomness and the angle cos _1 (p ■ p'/fc§), respec- 



tively. The coefficients (3j are expressed as (3^^+5(3^ 



where /3< wc) = (3^> = (3\^> = -(3^' = -2/3}" 
2/ g(wc) (T) = _^(2)( 1/2+ ( 4 ^TT)- 1 )/3 (r)/(7C(3)), and 
P (T) = 7C(3)A^(0)/(240tt 2 T 2 ). The strong coupling 
corrections 6(3j with impurity effects included were 
obtained, following Brinkman ct al. 6], in the spin- 
fluctuation model and, within the relaxation time ap- 
proximation, satisfies the relations 5f3^ — (5(3^ + 

5<5/^ sc) )/6, 5[3t c) = ^3 SC) + 55M SC \ and 5$ c) = 7<5/3< sc) , 
where 



,(wc) 



5f3[ sc) = -3.3 x 1Q- 3 f3 (T)tS ^{D (d \m)f, 



5(3t ] 



6D[ d \m)D { 2 d) (m) 



D {d \m) 



2(D{ d \m)f]/[J2( D i ) M) 2 ]> 

m 

[\rn\- 1 + (|m| + (2TTTT)- 1 )- 1 ' 



(3) 



\m\ + (AitTt)- 1 ) - t/>(l/2 + (4tttT)- 1 )}/2, D {d) (m) = 
V> (1) (1/ 2 + \m\ + (4tttT)~ 1 )/2, and the parameter S was 
defined in eq.(3.I0) of Ref.p. We find that, due pri- 
marily to the prefactor t = T/T c q in 6/3 j, the strong- 



2 



coupling corrections are weakened with increasing the im- 
purity strength (2ttTt) in qualitative agreement with 
the results in a static approximation 

In the relaxation time approximation used in obtain- 
ing eqs.Q) and the impurity-ladder vertex correc- 
tions are neglected by assuming 2ttTt » 1. Such vertex 
corrections are, at most, 0(l/(27rrT)) in magnitude and 
do not change the resulting phase diagram qualitatively, 
because they appear only as slight changes of a and /3q, 
i.e., in a manner independent of the pairing states. There 
are also additional terms 0, composed of a single im- 
purity line in /3m°^ leading to the change of coefficients, 
- -(-l) m A>eim P (m = 2, 4, or 5), where e imp 
is of the order (27tTt) -1 . The reason for our neglect of 
Eimp-terms will be explained later. 

To describe an inhomogenuity of An i due to the im- 
purity scatterings, an additional quadratic term 



n 2A = tY^KM^M)^ 



SjjSqj 2irN{0) 

3 2l£|+T- ] 



^g^KP^e(p + q/2,p' + q72) 
x £U(-p + q/2,-p' + q72) 



(4) 



usually neglected in the mean field (MF) analysis, needs 
to be evaluated, where Q £ (p, p') is the quasiparticle prop- 
agator defined prior to the random average and with a 
Matsubara frequency s. By expanding eq.Q) in pow- 
ers of u(r), its 0(u(r)) term identically vanishes, and its 
0([tt(r)] 2 ) term is evaluated as 



U 2A ~ T^^(q)^(q')^ 

<?,<?' e 



ir 2 k F N(0) 



2E F {2\e\+T^f 



A ' (2^^ Uk+qM - k - q '' (5) 

where we have assumed the k-integral to be dominated 
by large |k| values (see the 7's expression below). An ad- 
ditional quartic term arising from eq. (JHJ after the random 
average becomes 



K 



(n) 



10 



a, 6—1 qj 



x ^ 9 ,o(41fe)4!K?3))MW( 92 )^)(g 4 ) ) (6) 
where ^ 9 = 51+93 — 92 — Qi, a and b are replica indices , 



^ d= ^ E F T^42C(3) ^ (1) (2 ' 4.7V 



5 ^ 4 f.dDf 1 



1 



(7) 



and 7 = (tN(0)) 2 J dk 6(k F - k) 0{k) |w k | 4 /(27r 2 fc F ). 

We suppose that scattering events in an aerogel have 
local and strong anisotropics p), and that such random 
anisotropics leading to a pinning of the 1-vector in the 



ABM state be described as a strong fc-dependence of 
Uk| 2 . In fact, if |itk| 2 is a function only of fc 2 , and 
\q\, \q'\ <C k F , eq.© becomes an isotropic term like 
~ f r A* independent of the pairing state. Below, 

a strong random anisotropy such that |itk| 4 3> ( |uk| 2 ) 2 
is assumed. Further, since higher order terms in u(r) ne- 
glected above will also enhance 5 fid, the effective value 
of the parameter 7 may be large, e.g., of order unity, in 
3 He in an aerogel. Then, <5/3<j can become comparable 
with \8m sc ^\. Under these assumptions, the replicated 



GL functional ^ is given by w w = £ a (ft [A^] 



in) 



1(0)1 



^grl-^L"]]) + ^d' ij - H ^ S convenient to rewrite TL\ n> in the 



/(«) 



(») 



ABM state as H 



(n) 

d(i) 



(n) 
d(2) 



where 



(n) 
d(l) 



and 



(n) 
d(2) 



a, 6 



|A (a) (r)A w (r)| 2 , (8) 



|A^(r)A^(r)| 



x [ 1 - (l^ • l^) 5 



(9) 



For the BW and ESP states, 7Yd(2) does not appear, and 
the factor 3/5 in eq.(8) is replaced by 1/2. For a moment, 
we focus on the ABM pairing state. 

Below, the free energy will be examined in terms of the 
Gaussian variational method (GVM) 0, ^| which is of- 
ten used for random systems. By assuming a trial Gaus- 



sian functional 7i g 
GVM by 



(n) 



the free energy is approximated in 



F = F„ 



-(H 



(«) 

GL 



rL g 191 



where F g is the free energy for Hg 1 ^ divided by n, ( ) g 
is the ensemble average on TL g n \ and the n — > limit 
is taken at the end. It is reasonable to assume H. g n ^ to 



(10) 



take the form 7i a 



"^3,sym; where T^e^amp (^#,sym 



)is 



a Gaussian functional composed only of | A | -fluctuations 
(of symmetry variables such as the phase and /-vector). 
Then, the fact that Hgl appears in eq.^Hl only as its 
average permits us to use ea. (11UII after averaging over 1- 
orientations in terms including 9|A| in J2 a ^sxi^ul]- H 
can be accomplished with no knowledge on 7i ffiS ym in the 
present context with no strict LRO of the 1-orientation, 
i.e., 1=0. Consequently, Hgl can be replaced by the 



sum of the the amplitude part 7iim P 
where 



H. 



(n) 
amp 



H 



(n) 
d(l) 



= E /(«|A (a) | 2 + -f^(V|A^|) 2 +/3 A |A^| 4 ), 

(11) 



a=l • 



3 



and the symmetry variables' part Wsym 
where 



H 



(n) 



n 



O) 
d(2) 



L gl . 



2 2-j 

a— 1 



(|15|l and is the same as the T — > limit |nj of the cor- 
responding random-force model, which is obtained after 
the replacement, (l( a ) • 1») 2 -> 1 - (ft a ) - ft b )) 2 



/ |A| 2 (a^-)*(29j9 J -a^+V 2 a^ a ''.) (12) Ti^) m e q-03- Based on this fact, we will assume that, 



/3a = A245, /5ijfe = ft + ft + ft, and = \A\a^ t . The 
total superfluid free energy becomes the sum of -Fa mp and 

(n) (n) 

F sym , which are obtained below from H&mp and Wsym, 
respectively. 

Results on F amp are obtained by following the GVM 
analysis for the GL model of the Ising spin system [Icj . 
Since we focus on the GL region far below the critical re- 
gion, the relations (T/3 A ) 2 /[(iV(0)) 4 ^] « (T c -T)/T c0 < 
1 can be assumed together with (5ft <S ft. Then, follow- 
ing the analysis in Ref.^3] closely, F amp is expressed by 



F. 



amp 



N(0) \X P \ |A MF | 



V 



O(T/e ), (13) 



where V is the volume, the 0(T/£q) term is independent 
of the pairing symmetry, 



An 



3 5 / 2 T /3 A - 2<5ft/5 



N(0) (N(0We 



(14) 



corresponding near T c to (T — T c )/T c q, and |Amf| 2 = 
7V(0)|A p |/(2/3a)- Note that T c depends on the pairing 
state through the terms in eq. l|14|l except a/N(0). In 
fact, in clean limit with no (5ft, the relation /3a = /3b = 
fti2 + defining the bulk PCP easily follows un- 

der the condition that a state with higher T c is realized 
just below the resulting T C (P). Further, other fluctu- 
ation corrections of 0(T|A p |/^o) were dropped in F amp 
because, in the free energy difference between the ABM 
and BW states, the contribution from these corrections 
is of higher order in T / E-p compared with that from F sym 
(see ea. lfiHj) below). That is, F amp is well approximated 
by the condensation energy with a replaced by N(0)X p . 

Hereafter, ea. <|12|l will be treated in the London limit 
by replacing |A| there by |Amf|- Just as in F amp , the 
(5ft-independent fluctuation contribution to F sym will be 
neglected, and we focus below on the difference AF sym = 
-Fsym(<5/3d) — -Psym(O). However, a direct application of 
GVM to tiiym is not easy due to the 3D nature of L For 
this reason, the simpler model 



(n) 

sym(O) 



|A 



MF 



n 



{n) 
d(2) 



(15) 



will be first considered by assuming 1 to be a planar vec- 
tor (1 = x cos/3+y sinft) . For the model ea. ljTSj) . £±F sym /V 
in GVM was examined elsewhere to derive the quasi 
LRO and, when |AMF| 2 <5ft £,o/K t is small enough, is 
given by -3Tg c |A p |(5ftiV(0)/(40^ 2 X t /3A) < 0, where q c 
is a momentum cutoff of the order . This result is in- 
sensitive to the quasi LRO of 1-orientation in the model 



in order to evaluate the thermodynamic quantities, 7isym 
can also be replaced by the corresponding random-force 
model. Then, after a lengthy but straightforward calcu- 
lation, we obtain 



AF sym ((5ft is ) 
V 



-9 T\X P \ 
207t/3aCo 3 



<5ft. 



(16) 



For the present purpose, the sum of eas. (|13fl and (|16|l can 

be regarded as the total free energy in the ABM state. 

Here, the above results on the ABM state will be 
compared with those on other pairing states. The total 
free energy and A p in the BW ("robust" ESP) state are 
given by eq.@2J> an d ea. lCHfl if replacing /3a and 2(5ft/5 
there by /3 B (fti = ft + (ft 3 + 5ft 5 )/9) and <5ft/3, 
respectively, where ftj = ft + ft, while F amp and A p 
for the planar state are given simply by replacing /3a 
in them by /3p = P12 + ft45/2. Since it is easily seen 
that fti > /3p > /3a for any tT c q value we have exam- 
ined, higher T c (P)-values are obtained only for the ABM 
or BW states, and hence, the pairing states just below 
T c (P)-curve are the ABM or BW states. This result just 
below T c is unaffected by including the £i mp -terms in /3 m 
(see the paragraph prior to eq.|@J)). In fact, /3n is replaced 
then by /3n — fteimp (N = A, B, P, or R) commonly for 
the four pairing states. Further, the "robust" and 
planar states cannot occur even as an intermediate state 
between the ABM and BW states. This conclusion on 
the "robust" state with no AF sym is clear at this stage. 
We have also evaluated — AF sym for the planar state and 
found that it is at most T|A p |<5ft/(107r/3p£g), which can- 
not cancel out the cost in F amp relative to the ABM state. 

Examples of computed T C (P) and Tab(P) curves are 
shown in Fig.l. The T c (P)-curvc is defined from a larger 
value of T c for each of the ABM and BW states, while 
Tab(P) is obtained simply by comparing the free ener- 
gies of the two states with each other. In our computa- 
tion, we have assumed £i mp = 0. We have numerically 
verified that no visible change of Tab occurs by choos- 
ing £j mp = ±(2tttT)~ 1 . We have used Ef(P) data in 
Ref.[l3 and have assumed 5 = 300 T c0 /E F (P) together 
with a typical T c o (-P)-curve of bulk 3 He. The curve with 
(2ttt)^ 1 = 0.18 (mK) well explains why the A-like re- 
gion was not easily observed on warming |j| in spite of a 
significantly low PCP. Over a wide range of r _1 -values, 
with increasing r -1 , the width T C (P) — Tab(P) shrinks, 
while PCP tends to be lowered compared to the bulk 
PCP. The former feature is a reflection of a decrease of 
the strong-coupling correction to /3n, while the latter is 
due to an increase of <5ft. For low enough r _1 values 
(the (2-kt)- 1 = 0.13 curve in Fig.l), however, PCP is 
slightly higher than the bulk PCP. Such a nonmonotonic 
r-dependence of PCP is due to a competition between 
the two impurity-induced effects mentioned above. On 



4 





30 






la 


26 








24 






B 


22 








20 




18 




1.7 1.9 2.1 2.3 

Temperature [mk] 

FIG. 1: T C (P) and T AB (P) curves for (271-r)" 1 = (mK) 
(dotted curves), 0.13 (thick solid), and 0.18 (thin solid). Each 
PCP is denoted as a solid square on each T c (P)-curve. For 
all curves, 7 = 4.6 was commonly used. 



the other hand, a reduction of 7 (random anisotropy) 
under a fixed r monotonously shrinks T c — Tab and in- 
creases PCP. It is possible that a decrease of the poros- 
ity in aerogel corresponds to increases not only of t _1 
but also of 7. Then, the PCP in real aerogels may be 
monotonously lowered with decreasing the porosity. 

Finally, we comment on the ordering at long distances 
in the ABM state in aerogel in zero magnetic field by 
assuming a weak disorder at which no singular topolog- 
ical defects appear. Let us first consider the isotropized 
gradient energy H' glA + W^ (0) , where 



grA 



p s (vJ) 2 + A s (divl)^ + A b ((l- V)!) 1 



, (17) 



in place of Wsym for the ABM state. Here, A s , A/, > 



-IA 



MF 



"K t , and is the transverse component of ve- 



locity leading to nonsingular vortices 01 which might 
affect the ordering. In the limiting case with no first 
term, corresponding to the dipole-locked case in nonzero 
magnetic field where 1 is a planar vector (see ea.l|15|)). 
the quasi LRO of 1-orientation in the above model is al- 
ready known in the context of nematic glass 1 1 -* 3l| . To ex- 
amine effects of nonsingular vortices on this quasi LRO, 
the (vj) 2 term will be expressed using the Mermin-Ho 
relation as J r J r , fl(r) • fi(r')/(47r|r - r'|), where 
°c £ijkh{^lj x VZfc). The spatial nonlocality implies 
that this term is unrenormalized. Further, in the per- 
turbative renormalization analysis |13| under a fixed l 2 , 
this term quadratic in f2 is quartic in the "fast" variables 
like 1 1 x <51 1 . It implies that the (vj) 2 term is an irrel- 
evant perturbation to the T = fixed point controlling 
the l's quasi LRO. For the same reason, other vj-related 
terms, neglected in eq. G3, of niyl are also irrelevant. 
Thus, the l's quasi LRO is expected even in zero mag- 
netic field. Since the pair-field A^i is linear in d and 
the orbital "triad" including 1, the l's quasi LRO should 
result in the superfluid quasi LRO in any field fl5| . 

In conclusion, the free energy calculation including ef- 
fects of quenched disorder on the pair-field results in a 
phase diagram of 3 He in aerogel qualitatively consistent 
with experimental ones 0, and shows that the ABM 
pairing state describes the A-like state of 3 He in a erog el. 
Bearing the quasi LRO and a recent observation [l|| of 
an A1-A2 transition in mind, it is believed that the A-like 
state in 3 He in aerogel at weak disorder should be a glass 
phase with the ABM pairing and superfluid quasi LRO. 

We are grateful to LA. Fomin, K. Nagai, and T. 
Mizusaki for discussions. 
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